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ABSTRACT 

We present spectral analysis of five NuSTAR and Swift observations of GX 339-4 taken during a failed 
outburst in summer 2013. These observations cover Eddington luminosity fractions in the range ^0.9-6%. 
Throughout this outburst, GX 339-4 stayed in the hard state, and all five observations sho'w similar X-ray 
spectra 'with a hard power-la'w 'with a photon index near 1.6 and significant contribution from reflection. Using 
simple reflection models we find unrealistically high iron abundances. Allowing for different photon indices for 
the continuum incident on the reflector relative to the underlying observed continuum results in a statistically 
better fit and reduced iron abundances. With a photon index around 1.3, the input power-law on the reflector is 
significantly harder than that which is directly observed. We study the influence of different emissivity profiles 
and geometries and consistently And an improvement when using separate photon indices. The inferred inner 
accretion disk radius is strongly model dependent, but we do not And evidence for a truncation radius larger 
than 100 in any model. The data do not allow independent spin constraints but the results are consistent 'with 
the literature (i.e., a > 0). Our best-fit models indicate an inclination angle in the range 40-60°, consistent 
■with limits on the orbital inclination but higher than reported in the literature using standard reflection models. 

The iron line around 6.4 keV is clearly broadened, and we detect a superimposed narrow core as well. This 
core originates from a fluorescence region outside the influence of the strong gravity of the black hole and 'we 
discuss possible geometries. 

Subject headings: X-rays: individual (GX 339-4) — accretion, accretion disks — X-rays: binaries — stars: 
black holes 


1. INTRODUCTION 
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Black holes accreting via Roche-lobe overflow from a lo'w- 
mass companion star typically demonstrate a strongly tran¬ 
sient behavior: 'while they spend most of their time in qui¬ 
escence, during outbursts they can reach luminosities up to 
~ 10^^ erg s ', the Eddington li mit for a IOMq black hole 
( |Remillard & McClintock|2006| and references therein). The 
X-ray spectra can primarily be described by two components: 
a multi-temperature accretion disk and a hot electron gas 
corona. The corona shows up as a power-law component, 
originating from Compton up-scattering of soft disk seed pho¬ 
tons by coronal hot electrons. Outbursts typically start in a 
so-called low/hard state in 'which a power-law 'with a photon 
index T < 1.7 dominates the X-ray spectrum. The accretion 
disk contributes only 'weakly to the X-ray spectrum, with typ¬ 
ical disk temperatures kT ^ 0.2 keV. The disk is more evi¬ 
dent as a reflector of the hard X-ray continuum, resulting in 
a characteristic iron Kor line at 6.4 keV an d a Compton hump 
between 20^0 keV (|Ross & Fabian|20051l. In the later stages 
of the outburst, the source typically enters the high/soft state, 
in which the accretion disk gets hotter and becomes the domi¬ 
nant contributor to the X-ray spectrum, while at the same time 
the power-la'w becomes softer. 

In the high/soft state strong evidence exists that the accre¬ 
tion disk exte nds all the 'way to the inner most stable circular 
orbit (ISCO, Shakura & Sunya^|1973 i. The radius of the 
ISCO depends strongly on the black hole spin and 'with it the 
distortions of the fluorescent iron line due to orbital motions 
and general relativistic effects. By modeling these distortions, 
it is possible to estimate the spin of the black hole. An alter¬ 
native method to measure the spin is to use the continuum flux 
from the disk black body component, which is also a function 
of the ISCO dZhang, Cui & Chen|[T997) [Davis et al.||2005 
McClintock, Narayan & Steiner|2014[)' 
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In the low/hard state it has been postulated that the accretion 
disk is truncated and that the inner regions of the accretion 
flow are described by an optically thin, advection dominated 


outside of the low/hard state and obtain a - 0.93 + 0.05. |Reis| 


et al. (2008 1 , using XMM-Newton and RXTE data of both 
a soft and a hard state, confirm that value but require an in- 


accretion flow (ADAF, see, e.g., |Narayan & Yi 

[1995| [Esin,[ 

clination as low as i ^ 20° of the inner accretion disk. At 

McClintock & Narayan[[1997| Taam et al.[|2008 

1 . ITie trun- 

these very low inclinations fitting the disk continuum to mea- 


fore the fluorescent Fe Kor line from disk reflection should be 
narrow, and the blackbody radiation from the disk should be 
relatively cold. 

Measurements of the inner accretion disk radius at lumi¬ 
nosities > l%LEdd in the low/hard state typically are con¬ 
sistent with, or even require, an accretion disk extending 


2002 


all the way to the ISCO (see, e.g. Nowak, Wilms & Dove 
[Miller et aL]|2006[ [Petrucci et al.||2014[ [Miller et al. 


2014 i. ^eis, Fabian & Miller ( 2010| l extend these measure 


ments down to 0.05% Lsda using a sample of 8 sources and 
find that there is still evidence for an accre tion disk close to 
the ISCO. However, [Tomsick et al.j P009] l measure a trun¬ 
cated disk with an inner radius Ria > 175 for GX 339-4 at 
Lx ~ 0.14% Lndd, assuming an inclination of i - 30°. This is 
the lowest luminosity of GX 339-4 for which such a measure¬ 
ment is available and is the first time a clearly truncated disk 
is seen. The authors speculate that the accretion disk moves 
quickly outwards at luminosities <1% but it is not clear 
what mechanism triggers this change in the accretion disk ge¬ 
ometry. This idea of a recession below 1% Lndd in GX 339-4 


is supported by Allured et al. (20131 who measure inner radii 
around 10 for luminosities between 0.5-5% LEdd- 
GX 339-4 is the archetypical transient black hole binary. 
It has been studied in tensively since its discovery in 1973 
(Marker! et al. 1973|) and shows a high level of activity. 


with an average of one outburst of varying strength every 2 
years. During the hard-states it reaches fluxes in excess of 
150mCrab in the Swift/BAT 15-50keV energy band. This 
high flux and semi-regular outburst activity makes GX 339-4 
an ideal target to study the low/hard state in detail. 

The optical companion is not directly observable, and the 
only measurement of the orbital period has been obtained by 
measuring Doppler shifts of fluorescent lines c aused by the 
irradiation of the companion (Hynes et al.|2003|. These mea¬ 


surements give an orbital period P - 1.7 d and together with 
the upper lim it on the companion’s magnitude and assumed 
spectral type, Zdziarski et al. p004 1 estimate the black hole 


mass to be M — 10 Mq and the distance c/ = 8 kpc. These val¬ 
ues are^onsistentwithJhejrussJimitofM_^7Mocalculated 
by [Munoz-Darias, Casares & Martlnez-Paisj ( |2008| ) using a 
similar method and we adopt them throughout this paper. 

The inclination of the binary orbit is only weakly con¬ 
strained. As GX 339-4 is a non- eclipsing system, th e in- 
clination has to be lo wer than 60° (jCowley et al.|20()2 1 , and 


Zdziarski et al. (20041 est imate from the second ary mass func- 
tion a lower limit of 45°. jShidatsu et al. (|201 Ij) summarize all 


constraints and derive a best estimate of a; 50° 

The spin of the black hole in GX 339-4 is currently under 
debate in the literature. Estimates from reflection modeling 
and disk continuum models do not yet give consistent results 
regarding spin and inclination. This discrepancy is likely due 
to different underlying assump tions in the mod els (f or recent 
reviews on both methods, see |Reynolds||201^ and [McClin-j 
tock, Narayan & Steiner|20T4)l. However, both methods rule 


out a Schwa rzschild (i.e. non-spinning) black hole with high 
significance (Kolehmainen & Done|2010| Miller et al.|2008|l. 


Miller et al.j ( |2004||200}^| l use Suzaku and XMM-Newton 

data to measure the relativistic broadening of the Fe Ka line 


sure the spin gives consistent results, however, Kolehmainen| 
& Done (2010|) argue that this requires a strong misalignment 


between the orbit’s rotation axis and the spi n axis, as the or¬ 
bital inclinat ion is limited to be iorb >45°. [Kolehmainen &| 
|Done| ( |2()10| l therefore prefer a disk inclination > 45° which 
re sults in significantly lower spin (a < 0.9). 

[Plant, O’Brien & Fender] (|20T4| present XMM-Newton 
data taken during the same failed outburst in 2013 presented 
here (see below), and fit an inclination 33° + 3° while fixing 
the spin at a-0.9. They find consistent results between the 
disk continuum method and the reflection modeling assuming 
this inclination, and measure a small but significant trunca¬ 
tion of the disk, with an inner radius of 20-30 r^. They do 
not, however, address the misalignment between orbital and 
disk inclination. 

Therefore both the inner disk inclination and the spin value 
remain under discussion and need further investigation with 
sensitive X-ray instruments, especially covering the hard X- 
rays to get a go od measure of the con tinuum parameters. 

As shown by [Corbel et al.j ([2003[ and references therein), 
the radio flux of GX 339-4 is strongly correlated with the 
X-ray flux in the low/hard state, following the non-linear re¬ 
lationship f^adio This correlation has been connected 

to synchrotron radiation from the the jet also influencing the 


hard X-ray spectrum ([Markoff et al. [2003|l. [Corbel et al. 
( [2013[ ) have shown that the onset of in outburst is also ob- 
served in the radio spectrum, which switches from a negative 
spectral index a, where the radio spectrum is described by 
v", to a ^ 0.5 during the low/hard state. This correlation 
therefore provides information about the connection between 
accretion and ejection and between the jet and the corona. 

In August 2013, optical and X-ray monit oring detected 
the onset of a new outburst of GX 339-4 ([Buxton et al.| 


201 ^ 

the 


^ Pawar et al. 201 3j l. We triggered observations with 
Nuclear Spec troscopi c Telescope Array (N uSTAR, [Har- 


[rison et ar[2013| l. Swift ( [Gehrels et al. [20()4| ), and with the 
Australia Telescope Compact Array (ATCA). Overall we ob¬ 
tained five observations with NuSTAR, four during the rise 
of the outburst and one at the end of the outburst and Swift 
observations every other day, as well as three ATCA observa¬ 
tions. Figure [T] shows the lightcurve of the outburst, as seen 
with th e X-ray monitors Swif t/BAT ( Krimm et al.|2013 1 and 
MAXI ( [Matsuoka et al.|2009[ ). In a full outburst, the soft X- 
rays (as observed by MAXI) are expected to brighten as soon 
as the hard X-rays (as observed by BAT) decline, indicating 
the switch to the high/soft state. This outburst, however, did 
not follow that pattern and GX 339-4 never left the low/hard 
state resulting in a so-called failed outburst. Table gives 
a detailed observation log of the NuSTAR and simultaneous 
Swift observations. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Swift 

Swift monitored the complete outburst with ~ 1 ks snap¬ 
shots every other day. The evolution of the 2-10 keV count 
rate (scaled to m Crab units) as measur ed by the Swift X-ray 
Telescope (XRT, [Burrows et al.[[2005[) is shown in Figure [l[ 
For data reduction, we used HEASOFT v6.15.1 and the XCT 
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TABLE 1 

NuSTAR AND Swift observations and exposure times. 


No. 

ObsID NuSTAR 
(8000I0130XX) 

ObsID Swift 

MJD range 

exp. NuSTAR 
FPMA / B [ks] 

exp. XRT 
[ks] 

I 

02 

00032490015 

56515.907-56516.994 

42.25/ 42.26 

1.1 

II 

04 

00080180001 

56520.709-56521.914 

47.38/47.50 

1.9 

III 

06 

00080180002 and 
00032898009 

56528.525-56529.792 

43.78/43.94 

1.6 

IV 

08 

00032898013 and 
00032898015 

56538.414-56539.897 

61.94/62.29 

2.0 

V 

10 

00032988001 

56581.994-56584.299 

98.21 /98.22 

9.6 

Swift Illb 

N/A 

00032898008 

56527.788 

N/A 

0.97 

Swift IVb 

N/A 

00032898024 

56559.747 

N/A 

1.26 


200 



-20 0 20 40 60 80 100 


Time [d] since MJD 56501 

Fig. 1.— Light curves of SwiftfBAT (15-50 keV, orange crosses), MAXI (2- 
20keV, green circles), NuSTAR (3-79 keV, red diamonds), and Swift/XRT 
(2-10 keV, blue squares) of the 2013 outburst. All count-rates have been 
rescaled to mCrab fluxes in the respective energy band of the instrument. 
The right-hand y-axis gives the average measured NuSTAR count-rates of 
each observation. The downwards-pointing arrows at the top indicate the 
times of the ATCA observations. 


calibration released on 2013 March 13. For each ObsID that 
we used for the spectral analysis, we reprocessed the raw XRT 
data to produce new event files using xrtpipeline. Then, 
we used xselect to create source and background spectra. 
We included source counts from within 20 pixels (47") of the 
GX 339-4 position and background counts from an annulus 
centered on GX 339-4 with an inner radius of 90 pixels and 
an outer radius of 110 pixels. XRT was in Windowed Timing 
mode for the observations, and we scaled the background to 
account for the active detector area. For the response ma¬ 
trix, we used the file swxwt0to2s6_20010101v015.rmf and 
xrtmkarf to account for the effective area, including a cor¬ 
rection using the exposure map for each observation. The 
XRT spectra were rebinned to a S/N of 6 between 0.8-10keV 
using the Interactive Spe ctral Interpretation System (ISIS, 
Houck & Denicola||2000 1. All analysis was done with ISIS 


vl.6.2-17 in this paper unless otherwise noted and uncertain¬ 
ties are given at the 90% level. Observations Swift Illb and 
IVb listed in Table[T]were not used in the X-ray spectral anal¬ 
ysis as no simultaneous NuSTAR data are available. We use 
those observations only for comparison with the radio flux, 
as they occurred c lose st in time to the radio observations de¬ 
scribed in Section|23] 


2.2. NuSTAR 

NuSTAR consists of two independent grazing incidence 
telescopes which focus X-rays between 3-78 keV onto cor¬ 
responding focal planes consisting of cadmium zinc telluride 
(CZT) pixel detectors. NuSTAR, sensitive to X-ray energies 
from 3-79 keV, provides unprecedented sensitivity and high 


spectral resolution at energies above lOkeV, ideally suited to 
study the Compton reflection hump. The two focal planes are 
referred to as focal plane module (FPM) A and B. NuSTAR 
data were extracted using the standard NUSTARDAS v 1.3.1 
software. Source spectra were taken from a 120" radius re¬ 
gion centered on the J2000 coordinates. The background was 
extracted as far away from the source as possible, from a 135" 
radius region. This approach induces small systematic un¬ 
certainties in the background, as the background is known to 


change over the held of view and from chip to chip (Wik et al. 
2014|l. However, GX 339-4 is over five times brighter than 


the background even at the highest energies, so that these un¬ 
certainties are negligible. NuSTAR data were binned to a S/N 
of 36 in the relevant energy range of 4-78keV within ISIS. 
To reduce the spectral overlap which might be influenced by 
cross-calibration differences between Swift/XRT and NuS¬ 
TAR, we exclude NuSTAR data below 4keV. 

Timing analysis of Nu STAR observations I-IV is presented 
by [Bachetti et al.| ( |2015[ ). The power spectral density (PSD) 
between 0.001-200 Hz is consistent with typical hard-state 
PSDs an d well described by three Lorentzian components 
(see, e.g. |Belloni et al.||2005||. We also calculated the PSD 
for observation V and And consistent results. 


2.3. ATCA 

We obtained quasi-simultaneous radio observations with 
ATCA, as a radio jet is expected to be launched during the 
low/hard state of GX 339-4 ( [Corbel et al.1[2000l l. ATCA’s 
synthesis telescope is located in Narrabri, New South Wales, 
Australia, and consists of six 22 m antennas in an east-west ar¬ 
ray, using linearly orthogonal polarized feeds which allow the 
recording of full Stokes parameters. The observations were 
conducted simultaneously at 5.5 and 9 GHz on MJD 56516.4, 
56527.46, and 56559.29, i.e. simultaneous to observation I, 
Swift Illb, and Swift IVb (see Fig. [H, using the u pgraded 
Compact Array Broadband Backend (CABB) system (]Wilson| 
et al.|20lT]l. The first observat ion was conducted and reported 
by |Miller-Jones et al.| ( |2013[ ). The array was in a compact 
configuration (H214 or H168) during this period of ATCA 
observations. 

The amplitude and band-pass calibrator was 
PKS 1934-638, and the antennas’ gain and phase cali¬ 
bration, as well as the polarization leakage, were usually 
derived from regular observations of the nearby calibrator 
PMN 1646-50. The editing, calibration, Fourier transfor¬ 
mation with multifrequency algorithms, deconvolution, and 
image a nalysis were performed using the MIRIAD software 
package (jSault & Killeen|1998|l. 

We show the observed 9 GHz fluxes as a function of the 3- 
9 keV X-ray flux as measured with Swfft/XRT in Figure 
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Fig. 2.— ATCA 9 GHz flux as function of Swiff/XRT 3-9 keV flux. Our 
2013 data are shown with red diamonds, data from failed outbursts in 2008 
and 2009 presented by|Corbel et al.|j2013) as black circles and black squares, 
respectively. The superimposed black line is the be st-fit correlation from all 
archival data (containing full and failed outbursts, [Corbel et al.|2013t , the 
orange line is the best-fit to the observations during tailed outbursts. Note 
that the error-bars on the 2013 data are smaller than the symbol size. 


and compare them with archival data presented by Corbel 
|et al.| ( |2013| l. The data points fall below the measured cor¬ 
relation from all outbursts given by these authors and agree 
better when taking only the data of other failed outbursts in 
2008 and 2009 into account. This behavior seems to indicate 
that the jet power is somewhat reduced in failed outbursts and 
hints at a different accretion geometry. 

3. SPECTRAL MODELING 
3.1. Basic Fits 

As shown in Fig. all five observations show very similar 
spectral shapes in NuSTAR, with clear evidence of reflection 
features. To highlight these features, we fit each NuSTAR ob¬ 
servation with a simple absorbed power-law, using only data 
between 4-6 keV, 8-10 keV, and 50-78 keV, i.e., ignoring the 
energy ranges where the strongest contribution from reflection 
features is expected. The residuals shown in Fig. [^clearly in¬ 
dicate a strong Fe Ka line and Compton hump. The shape 
of the iron lines appears to be constant over all observations, 
while we see indications that the Compton hump is more sig¬ 
nificant in the high flux data. 

For the remainder of the paper we fit all Swift/XRT data 
between 0.8-10keV and all NuSTAR data between 4-78 keV 
simultaneously, unless otherwise noted. To model the reflec- 


tion component we add the reflionx model (Ross & Fabian 
[2005)l, with its high-energy cutoff energy set to 500 keV. 
While this model was originally calculated for AGN spec¬ 
tra, it has been applied successfully in the literature to de¬ 
scribe black -hole binaries (see, e.g., the discussion in [Fabian] 
jet al.||2012| l. This model self-consistently describes the ef¬ 
fects of reflection of an input power-law spectrum off an op¬ 
tically thick accretion disk, including the Compton hump and 
the fluorescent iron lines. For consistency we initially tie the 
photon index of the reflionx model to that of the primary 
power law (model Ml). To allow for relativistic effects close 
to the black hole we fold the model components th rough the 
relconv convolution kernel (jPauser et al.||2010 1 . We fix 


the outer radius of the accretion disk to Rout = 400 r^, as at 
these distances relativistic influences are negligible and the 
reflected flux is expected to be low. 

In XSPEC notation the model is represented as constant 


Fig. 3. — (a) Unfolded spectra of all five observations. Only NuS- 
TAi?/FPMA data are shown for claiity. Observation I is shown in red, ob¬ 
servation II in blue, observation III in magenta, observation IV in orange and 
observation V in green, (b) Ratio residuals to a simple power-law fit. See text 
for details. 

* tbabs * [powerlaw + relconv (reflionx)]. Here, 
constant takes into account small flux differences in the cal¬ 
ibration of the different instruments. We quote all fluxes rela¬ 
tive to NuSTAR/FPMA. 

While the continuum and reflection parameters are ex¬ 
pected to change over the outburst, some parameters can be 
assumed to stay constant. In our model, these include the 
Galactic absorption towards the source (A^h). the black hole 
spin (a), the inclination of the accretion disk (/) and the iron 
abundance in the disk (Fe/solar), expressed in units of solar 
iron abundance. To obtain the best statistics for these parame¬ 
ters, we fit all five observations simultaneously, requiring that 
these parameters are the same for all observations. 

The absorpti on is described using an updated version of the 
tbabs model ( Wilms, Allen & McCray |2000| l, with the cor- 
responding abundances and cross-sections fromjVerner et al.| 
((1996)). The fitted absorption column is only marginally 
higher than the expected Galactic a bsorption column of 
3.74 X 10^' cm“^ ( Kalberla et al.|2005 1 and is in good agree¬ 
ment with values used in literature for this sour ce (e.g., Tom- 
jsick et al.|2009(|Plant, O’Brien & Fender|2014 1 . We therefore 
do not add an additional absorption column intrinsic to the 
source. _ 

As a first approach (model Ml-q3, Table |AT|), we assume 
a standard Shakura-Sunyaev accretion disk jshakura & Sun- 
yaev 1973|l, and we use an emissivity index (q) of 3 ( [Reynolds 
& Begelrnan[[1997[). The data quality does not allow us to 
constrain the spin of the black hole at the same time as the 
inner radius of the accreti on disk. We therefor e fix the spin to 
a = 0.92, as measured by [Miller et al.|([2008]). As described 

[flu 


in Sect, 
niflcant 


4.2| the choice of a does not influence our results sig- 

y- 

Allowing the emissivity index q to vary individually for 
each observation (model Ml-qv, Table [A2] l improves the fit 
1698 dof and results in a steeper emissivity 
~ 5 for all observations but observation V, where it is fitted 
to q '5 = 1.7 + 0.5. All other parameters do not change signifi¬ 
cantly. 

The relconv model also allows us to use the “lamppost” 
geometry (Ml-LP, Table [A3|l, where the corona is assumed to 


be a point source on the spin axis above the black hole (Mini 


utti & Fabian 2004[). The main free parameter in this geometry 
is the height of the corona above the accretion disk plane and 
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the illumination and emissivity index of the accretion disk is 
self-consiste ntly calculated takin g effects of general relativity 
into account (|Dauser et al.|201() i. This model also provides a 


small improvement above the Ml-q3 model - 1-37 for 
1698 dof) and the coronal height is fitted to be < 4.2 for 
all observations. All other parameters do not change signifi¬ 
cantly. 

Including a soft black-body component in any of the models 
does not improve the fit quality significantly nor does allow¬ 
ing the cutoff energy to vary. The 90% lower limits on the cut¬ 
off energy are, in all observations, > 370 keV. Using models 
not accounting for a Compton hump, [Miyakawa et al. ( |2008 1 
and |Motta, Belloni & Homan | ( |2009[ ) hnd somewhat lower 
cutoff energies at similar luminosities (0.9-5.6x10^^ erg s“' 
between 2- 200 keV). However, wh en including a reflection 
component, [Miyakawa et al^ ( |2008| l do not find evidence for a 
cutoff at E < 500 keV at these lurninosities . This is in agree¬ 
ment with results presented by [Plant et al.|(|2014|l. The con¬ 
tinuum in our data is therefore \vell described by a power-law, 
consistent with previous studies. 

In all models (Ml-q3, Ml-qv, and Ml-LP) we find an 
unphysically high iron abundance (5.00 ^q } 2 , 6.5 + 0.4, and 
5.27^g2g solar, respectively). In the literat ure, abundances of 


1 or 2 times solar are typically assumed (Miller et al. 2008[ 


Reis et al.|2008[[Tamura et al. [2012[[Allured et al.|2013 [Plant| 

et al.|2014[ l. Forcing a lower iron abundance results in clearly 
worse hts, with strong residuals around the Fe Kcr line en¬ 


ergy. Using the xillver reflection model (Garcfa & Kall- 
manj2010| instead of reflionx does not improve the ht and 


requires an equally high or higher iron abundance. 



Time [MJD] 

Fig. 4.— Results of the spectral fit for models M2-q3 and M2-LR Results 
of model M2-q3, i.e. with a fixed emissivity q = 3, are shown as black 
diamonds, while results for model M2-LP, i.e. using the lamppost geometry, 
are shown as red circles. The latter are shifted in x-direction for clarity, (a) 
photon index of the primary continuum, (b) ionization parameter, (c) input 
photon index for the reflection component, (d) inner disk radius in rg. 


3.2. More sophisticated modeling 
3.2.1. Reflector sees a different continuum 

The large iron abundance can be lowered and the statisti¬ 
cally unacceptable values improved by using more com¬ 
plex models, going beyond the traditional power-law plus 
reflection model. The standard geometries assume that the 
corona is a point-source and uniform, but we have to expect 
that the geometry in reality is more complex. The simplest 
approach to describe a corona which is physically extended 
with a non-uniform temperature profile is to allow different 
photon indices for the observed power law continuum and the 
input spectrum to the reflionx reflection (model M2). 

When fixing the emissivity index to q - 3 (M2-q3), this 
model improves the x^ value significantly by Ax^ = 268 for 
five additional parameters over Ml-q3. We show the best- 
fit parameters for this model in Table 1^ We also apply the 
lamppost geometry (M2-LP) for which we give the best-fit 
values in Tablej^ Both models give mostly similar results and 
a comparable quality of the fit. The evolution of the spectral 
parameters with time for both models is shown in Figure]^ 

We show the residuals for the M2-q3 model for all five ob¬ 
servations in Figure As can be seen, there is good agree¬ 
ment between all three instruments, though the NuSTAR data 
clearly dominate the statistics. Small systematic residuals at 
the lowest NuSTAR energy end can be attributed to calibra¬ 
tion uncertainties. All five observations show very similar 
residuals and have a comparable quality of the fit. 

The inner radius is highly unconstrained, with best fit val¬ 
ues around 100 rg for M2-q3 but close to the ISCO for M2- 
LR Both models are marginally consistent with values around 
50 rg. Due to the large uncertainties and geometry dependence 
of the values, it is not clear if a truncation of the accretion 


disk is present or not. In the lamppost geometry the large un¬ 
certainties can be understood as being due to the large coro¬ 
nal height, which results in an emissivity index clearly be¬ 
low 3 between ai3-100rg ([Vaughan et al.|2004 Dauser et al.] 
[2010j l. Therefore, the inner parts of the disk contribute less 
to the reflection model and the inner radius cannot be well 
constrained. 

We find very similar coronal heights around 150 rg for each 
observation with large uncertainties. The height is completely 
unconstrained in observation V in the allowed range between 
3-300 r^. We therefore conclude that we see no indication 
for an evolution of the corona over the outburst and that the 
corona seems to be located relatively far away from the black 
hole. 

In the M2-q3 model the inclination is fitted to be i — 48°;^'^ 
(49° in M2-LP). These values are in very good agreement 
with the limits on the orbital inclination so that no misalign¬ 
ment between the accretion disk and the orbital plane is nec¬ 
essary. 

We also obtain similar results when allowing the emissiv¬ 
ity index to vary (M2-qv). The best-fit emissivity indices are 
relatively flat, around 1.6-1.8 and the models indicate a non- 
truncated accretion disk. The other parameters do not change 
significantly. The fit statistics improve by Ax^ - 46 for five 
additional degrees of freedom over the M2-q3 model, which 
correspo nds t o an F-test false alarm probability of 6 x 10“^, 
see Table lA4l 

The reflection fraction R given in Tables [^ and [^ is the ra¬ 
tio of the flux of the reflection componentbetween 0.001- 
lOOOkeV to the flux of the primary unabsorbed power-law 
component between 0.1-1000 keV. While this implies extrap¬ 
olation of the model far outside the fitting range, it captures 
the whole energy range used in the calculation of reflionx. 
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TABLE 2 

Best-ht parameters for emissivity index q = 3, spin a = 0.92 and allowing for two different photon indices (M2-q3). 


Parameter 

I 

II 

Ill 

IV 

V 

Ah [1022cm-2] 

0.868 ± 0.020 

_ 

_ 

_ 

_ 

Fe/solar 

1 7T+009 
^•'■^-0.08 

— 

— 

— 

— 

i [deg] 

48+‘2 

A rr^^n+O.OOU)- 

— 

— 

— 

— 


A power 
Trefl 

Rin [rg] 
CCxRT 


1.585 ±0.009 
1 


(0.97+O- 


-0.05 

|3)xl0-5 

2.21+«;22) X 102 

1.7^1 ^)x 102 
0.945 ± 0.028 


1.312+ 


1.594 ±0.007 
+0.025 
-0.033 

(i.67+;J;|2)x10-3 

o.iol^,n2 
0.08 j 


0 . 10 , 

228+'>-1'>' 


') X 102 


0.9+1;^) X 102 

0.915 ±0.015 


1.616 ±0.006 
1 ooT+0.017 

i.o:)j_ 0 .oi 6 
(2.25 ± 0.14) X 10- 
2.65+“ ;^) X 102 
1.3+* *) X 102 
1.129 ±0.014 


1.643 ±0.005 
7+0.016 
-0.015 

;.24+“|S)x 10-3 
2.52+“+“ X 102 


1.357 
(3, 


0.65+“ 32) X 102 

1.104 ±0.012 


1.608 ±0.007 
1.34 ±0.04 
(4.2 ± 0.4) X 10-3 
2.34!“;|2)x102 

2.3^“ *) X 102 

1.052 ±0.014 


AT 

0.30 ± 0.06 

0.282 ± 0.029 

0.284 ±0.018 

0.285 ±0.016 

0.26 ± 0.04 

R 

0.64 

0.71 

0.78 

0.83 

0.52 

%iedd’’ 

1.65 

2.74 

4.35 

5.63 

0.93 

^2/d.O.f. 

2124/1698 






1.25 





^ in phs“^ cm“^ 






^ Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of lOM© 




TABLE 3 




Same as Table[^but for the lamppost geometry 

(M2-LP). 


Parameter 

I 

II 

III 

IV 

V 

Ah [I 022 cm- 2 ] 

0.85lt“2“ 

— 

— 

— 

— 

Fe/solar 

1 7 A+O.O 6 
^•'^-0.08 

— 

— 

— 

— 

i [deg] 

49+2 

— 

— 

— 

— 

'^conl^ 

0.0668f“;““i“ 

A 1 1 ao+0.0013 

noi? 

A 1 OAQ+0.U016 

U.18U9_qqq20 

^•^'^^^-0 0017 

A A^i O+U.UU07 

r power 

1 coc+0.008 

i.JOJ_QQ06 

1 T 1 +0.05 
-0.06 


i.uij -0 006 

i.ohuz_qqqj 3 

1.605 ±0.006 

Trefl 

1 714+0.025 
i.2)i‘+-0 028 

1.336 ±0.015 

1.363)“;“!“ 

1.35 ±0.04 


(i.02)““8 X 10-3 

(1.71 ±0.10) X 10" 

-3 (2.36)“ |“) X 10-3 

(3.32)“;^ X 10-3 

(4.2 ± 0.4) X 10-3 


( 2 .IF“;W)x 102 

(2.23)“;“«)x102 

(2.54)“;‘3)x102 

(2.45)“;““) X 10+ 

(2.29)“;‘“)x102 

Rin [rg]*^ 

< 66 

< 74 

<59 

< 26 

< 300 

//(corona) [r^] 

(i.5)‘;3)xI02 

(1.4)+“) X 102 

(l.9)+“)x 10+ 

(l05)“|2)x102 

( 3 . 00 )“““) X 10+ 

CCxRT 

0.946)““23 

0.914)““12 

1 IOQ+9.012 

1-I^o noi5 

1.103 ±0.010 


AT 

0.28 ± 0.05 

0.278 ± 0.027 

0.278 ±0.016 

0.278 ±0.015 

0.26 ± 0.04 

R 

0.63 

0.70 

0.77 

0.81 

0.51 

LJLedd X 100“ 

1.63 

2.73 

4.32 

5.60 

0.93 

V/d.o.f. 

2103/1693 





xid 

1.24 






‘ in ph s ^ cm ^ 

’ Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of IOMq 
' The lower limit of Rm is the ISCO at 2.2 rg for an assumed spin of a = 0.92. 


The reflection fraction is correlated with the X-ray luminosity, 
as already indicated by the variable Compton hump strength 
shown in Figure]^ As the inner accretion disk radius seems to 
be constant, this indicates a change in coronal geometry where 
at higher flux more of the coronal emission is intercepted and 
reprocessed by the accretion disk. 

The reflection fractions given in the tables are no t directly 
comparable to the one given by |Plant et al.| (|20T4]l, as these 
authors use only the 4-lOkeV energy band to calculate it for 
the outbursts between 2002-2008. In this energy band we 
measure a refl ection fraction aro und 0.05, about a factor of 
2-3 lower than Plant et aL] ( |2014| l. 

The photon index incident to the reflection component is 
very hard (Frefl ~ 1.3), indicating a strongly photo n starved 
Comptonization region (|Haardt & Maraschi||1991|l. Its dif¬ 
ference to the photon index of the primary power law, AF, is 
relatively constant over the outburst in all models; AF ^ 0.3. 


If this difference is due to the geometry of the corona, the 
corona seems to stay stable, even though the reflection frac¬ 
tion changes as a function of flux. We note, however, that 
the lamppost intrinsically assumes a point-like corona so that 
the different photon-indices we find here cannot directly be 
geometrically interpreted. 

If we want to model an extended corona along the spin 
axis, perhaps corresponding to the extended base of a jet 
( jMarkoff, Nowak & Wilms.2005| l, this can be done more self- 
consistently in XSPEC or ISIS by adding a second relativis- 
tically smeared reflection component, with a coronal height 
larger than the first component and with its photon index tied 
to the observed primary continuum (M2b-LP). That is, we re¬ 
quire that the observed continuum is also reflected and add as 
additional free parameters the flux of its reflection component 
and the second lamppost height. The incident spectrum to 
the first reflector is understood to originate much closer to the 
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0.8 1 ^_^^. . .^^^ 

1 2 5 10 20 50 

Energy [keV] 


Fig. 5. — (a) Data and best-fit model for M2-q3 for observation IV, show¬ 
ing the reflection component in magenta and the power-law components in 
orange. Swift/XRT data are shown in green, MiSTAR/FPMA in red, and 
FPMB in blue, (b)-(f) Residuals in terms of model to data ratio of the best-fit 
model for observation I-V, respectively. Data were rebinned for clarity. See 
Table|^for the model parameters. 

black hole such that it is invisible to the observer. This model 
results in a good ht, with =1.214 which corresponds to 
an improvement Ax^ - 54 for six additional degrees of free¬ 
dom over the M2-LP model. While this is only a marginal 
improvement, the added self-consistency is important for the 
physical interpretation of the model. 

Table 0] shows the best fit values of this model. As in the 
other M2 models, the inner accretion disk radius is consistent 
with being at the ISCO. The inner reflection height. Hi is fit 
to the minimal allowed value, 2.2 r^, while the outer reflection 
height, H 2 , pegs at the upper limit, 300 r^. Thus the model 
describes a very strongly elongated corona, and one would 
expect that intermediate heights also contribute to both the re¬ 
flection and the observed primary continuum. Such a model, 
however, is not uniquely definable in ISIS or XSPEC as all re¬ 
flectors and continua would be degenerate. We therefore take 
the two reflectors presented here as the best approximation to 
an elongated corona with varying power-law emission. 

3.2.2. Secondary continuum component 

In continuation of the idea that two different hard power- 
laws are present in the system, perhaps from gradients in the 
coronal temperature, we apply a model which consists of two 
power-law components, with only one being the input to the 
reflection component (M3). We tie the photon index of one 
power-law to the photon index of the reflection component. 
With respect to the models presented in the previous section 
(M2), we have another free parameter, the normalization of 
the second power-law component. 

We again fit three different models: the first with the emis- 



Fig. 6 .— Same as Figure|^but for model M3-q3. The power-law incident to 
the reflection model is shown in orange, the second power-law in blue. Resid¬ 
uals are shown in units of cr. The best-fit parameters are given in Table|^ 


sivity fixed at q' = 3 (M3-q3), the second allowing the emis- 
sivity index to vary (M3-qv), and finally using the lamppost 
geometry (M3-LPX Figure shows the residuals for M3-q3, 
separately for each observation for clarity. Table gives the 
best-fit values for M3-q3. We give two values for the reflec¬ 
tion fraction: R\ is calculated using only the first continuum, 
which has the photon index used for the reflection component. 
R \+2 is calculated using the total observed flux from both con¬ 
tinuum models. 

All three models result in very comparable qualities of fit 
with TL/dof= 1.22/1693 (M3-q3), 1.21/1688 (M3-qv), and 
1.22/1686 (M3-LP). However, the relative strengths of the 
two continuum components depend on the assumed geometry. 
The lamppost model, in particular, gives ver y dif fer ent r esults 
(all model parameters can be found in Tables A5 and A 61 . Ad¬ 
ditionally, in observation V, the first power-law norrnalization 
in all three models is cons istent with 0, reducing this model to 
the M2 described in Sect. |3. 2.1 In the M3-q3 model the sec¬ 


ond, non-reflected power-law contributes between 30-50% of 
the broad-band flux in observations I-IV. 

The best fit iron abundance is found to be between 2-2.5 
the solar value, depending on the geometry. The inner radius 
of the accretion disk is again consistent with being close to 
the ISCO, especially in the model with a free emissivity in¬ 
dex where the best fit values are below 20 Cg. The differences 
between the two photon-indices are typically larger than in 
the previous model, with values around AT = 0.5. Overall, 
this model confirms the previous results that we can obtain a 
reasonable iron abundance by using a more complicated con¬ 
tinuum model and that the accretion disk is not significantly 
truncated. However, through the introduction of another free 
parameter, the parameters are less well constrained and the fit 
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TABLE 4 

Same as TableI^but eor the lamppost geometry and two reelection components whose input spectra originate erom dieferent 

HEIGHTS (M2b-LP) 


Paiameter 


III 


IV 


Ah [10^2 

Fe/solar 
Al (bottom) [rg] 
A2(top) [rg] 
i [deg] 
A 

^conl 

Tpower 

Treli 


CCxRT 


0.818 ±0.016 
2 

^■■^“-0.27 
9 9 + 1.0 
^•^-0.0 
(3.00+®®°^ 


X 10^ 


-1.44, 

39.2+2;’ 
-QQZi^+truon” 

U.U04J_.pQj^ 

i.juo_oo09 

1 10 +®*^ 
‘■++'- 0.10 


0.1075 ±0.0014 

1 ^on+0.008 
i.jou_goo7 
I 1 q+0.08 
^•^^-0.10 


0 1755+®-®®*** 

+'■+ '■'+-0.0017 
1.602 ±0.006 
1 1 O+0.09 

1.1»_0 17 


0 2461+®-®®’® 

U.4±D1 QQQJg 

1.630 ±0.005 

1 16+®-®* 
'■"'-0.07 


^Lrefl" 1 

J0.67+0;[6) X 10-6 1 

''l 04+0.16'' 
-0.14/ 

1X 10-6 1 

''1 09+0.25'' 

t^-^^-0.18/ 

1 X 10-6 

(2.80/0;«)x10-6 ( 

■^2,refl^ 1 

4.026+® ®*®) X 10-6 1 

;o^5oti4 

1X 10-6 1 

u.oz_0 15 , 

1 X 10-6 

(1.05 ± 0.20) X 10-6 ( 

^ 1 

J2.13!«;‘6)xl02 1 

0 90+0.10' 

1 X 10’ 1 

7 cQ+0.20' 

r-^^-o.i6> 

0 

X 

(2^429!«®^*)x102 ( 

[rg] 6.3/*®^’ 5.2+[;6 

4 3+* -^ 


3.2 ±0.6 ( 


0 941 +0 028 
"■++'-0 027 


0.914 ±0.015 


1.131 ±0.014 


1.104 ±0.012 


0.0407 ± 0.0007 
1.590 ±0.008 
1 •! 1 +0.10 
-0.31 


3.54+ 


+0;*®) X 10-6 
X 10-6 
X 10’ 

X 10’ 
1.053 ±0.014 


2.26!«;|^) 


AT 

0.47 ±0.12 

0.39 ±0.10 

0.43 ±0.13 

0.47 ± 0.07 

0.28 ± 0.20 

R 

0.63 

0.69 

0.77 

0.84 

0.47 


1.65 

2.73 

4.34 

5.62 

0.94 

X^ld.o.f. 

Xli 

2053/1691 

1.21 






‘ in ph s * cm ^ 

’ Luminosity calculated between 0.1-300 keV assuming a distance of 8kpc and a black hole mass of 10 M© 


TABLE 5 

Same as Table[2but for 2 power-law continua (M3-q3). 


Parameter 

I 

II 


III 

IV 


V 

Ah [10 ’’cm-2] 0.904+o^022 

Fe/solar 2.07+® 

I [deg] 31+6 

= 



E 



'^cont,!^ ^ 

rcont,! and refl 

^refl 1 

Rm [rg] 1 
'4cont,2^ ^ 

rcont,2 
CCxRT < 

1 034+®®®® 

■"+ -Q.Q14 

1.426+|i;i 

;0.59;®;«®)x10-6 

*2.51+®;’6 X 102 

;i.55/®;®«)xl0’ 

3 rj97+0 010 

j.uj /-Q po7 

1 Q 1+0.24 
^•^^-0.18 

3.93 ± 0.04 

0 ( 

1 420+^^^* 

-0.036 

(l.ll+«|6)xl0-6 1 

(2.49+« 22) X 102 1 

44^2® 1 

0.066;;®^® ( 

1 Q9+0.17 
-0.14 

0.900 ±0.018 

3 075+® ®*6 

4 ;;8:8?§ 

-0.019 

;1.52+®;2®)x10-6 

;3.03!0;66 X 102 

;0.52!0;«)x102 

\ 1 1Q+0.016 
-'••^^^-0,010 

1 89+013 

^•°^-o.io 

1.120 ±0.017 

~Q~Q^^TTnm 

1 47JO+0.015 

1 .-+JJ_0.020 
(2.3 ± 0.4) X 10' 
(2.85+“6) X 10* 
26/)6 

0 171+0.027 
'^-Q.pl6 

1 87+01? 

-0.09 

1.100 ±0.013 

( 

1 

^ 1 

( 

TQQQ^yhrnTTTTT 

/ _o 00027 

1 oQ+0.05 
^•'^^-0.04 
)3.6/®^^) X 10-6 

*2.53+«:22)x102 

,0.7!®;:) X 102 

3 041 7+0 0009 

1.603lO;Jg 

1.068ii 

AF 0.48 ±0.21 

0.50 ±0.16 

0.47 ±0.12 

0.44 ±0.10 

0.22 ± 0.08 

Ri 0.46 

0.52 

0.62 

0.74 

29.87 

Ri+2 0.31 

0.33 

0.36 

0.39 

0.49 

%Ledd6 2.16 

3.66 

5.88 

7.71 

0.94 


A^/d.o.f. 2072/1693 

1^22 


" in ph s ' cm ’ 

6 Luminosity calculated between 0.1-300 keV assuming a distanee of 8 kpe and a blaek hole mass of IOMq 


quality is not significantly improved. 

3.3. Iron line complex 

In all previous fits we modeled the Fe Ka line self- 
consistently as arising from reflection off an accretion disk 
with constant ionization, as described by the reflionx 
model. However, close inspection of the residuals reveals that 
none of these models captures the shape of the line perfectly. 
We therefore add an ad-hoc narrow Gaussian component to 
the M2 models, with its energy hxed at 6.4keV, the energy 
of neutral Fe Kcr. This line could, for example, be produced 
further out in the disk, where we encounter near neutral iron 
and the influence of the relativistic effects of the black hole 


are negligible. With this addition, we hnd signihcantly better 
hts for all versions of the M2 model, i.e., when the photon 
indices between the continuum and reflector are i ndependent. 

We note that [Plant, O’Brien & Fenderj ( |2014| l do not see 
evidence for such an additional narrow component. How¬ 
ever, their data lack the crucial coverage above lOkeV and 
therefore the reflection model is entirely driven by the iron 
line shape. When applying a model similar to their best-ht 
model using the lamppost geometry to our observation IV, 
and limiting the htting range to below lOkeV, we obtain very 
similar results, most notably an ionization around ^ = 1000 
and no residuals in the iron line band. However, this model 
clearly underpredicts the flux in the hard X-ray band, espe- 
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TABLE 6 



Same as Table[^but adding an additional narrow Gaussian line at 

6.4 keV (M2-q3-Ee). 


Parameter 

I 

II 

III 

IV 

V 


n 

1 CQ-1-0.10 

^•'’°-0.09 

— 

— 

— 

— 

Fe/solar 

— 


— 

— 

( [deg] 

59/*^ 

— 

— 

— 

— 


0.0667 ± 0.0010 

O.1114 ±0.0014 

0.1823 ±0.0018 

0.2556 ± 0.0023 

0.0418 ± 0.0006 

r cont 

i.joj_q009 

1.594 ±0.007 

i.ui / /_0 0023 

1.645 ±0.005 

1 603 

-0.004 

r,efl 

1.24 ±0.06 

1.27 ±0.04 

l.i»U_oo27 

1 tiq+0.017 

1.25 ±0.06 


(l.08/‘“/)x 10-5 

(1.75 ± 0.12) X 10-5 

(2.40 ± 0.16) X 10-5 

(3.45/«;|2 X 10-5 

(4.4 ± 0.5) X 10-® 


( 2 . 10 /°‘^jx io2 

(2.24/0;®) X 102 

(2.58/0;i^)xl02 

(2.44/«;>«)xI02 

(2.31/«;|5)x102 

Rin Fj] 

TQ+34 

(0.72/0;«)xI02 

51/50 

40+20 

-13 

45/« 

ApeKti^ 

(0.77/[J-2^)xI0-^ 

(0.7 ± 0.4) X lO-'^ 

(1.7/0 ^) X10-" 

(2.1 ±0.5) X lO-'* 

(0.51 ±0.13) X 10-"^ 

CCxRT 

0.949 ± 0.028 

0.912 ±0.015 

1.124 ±0.014 

1.099 ±0.012 

1.052 ±0.014 

AT 

0.34 ± 0.06 

0.33 ± 0.04 

0.338 ± 0.027 

0.327 ± 0.017 

0.36 ± 0.06 

R 

0.65 

0.73 

0.80 

0.85 

0.53 

Eqw (FeKa) [eV] 

22/2 

12 ±7 

19 ±5 

17 ±4 

24+« 

%Z.cdd 

1.64 

2.74 

4.36 

5.62 

0.94 


2009/1693 






1.19 






^ in ph s * cm ^ 

^ Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpe and a black hole mass of 10 M© 


dally around the Compton hump. 

As an example of a model with an additional Gaussian com¬ 
ponent, we give the best values for a fixed emissivity index of 
^ = 3 in Table (M2-q3-Fe). The narrow core is fitted to 
an equivalent width around 20 eV in all observations and does 
not show a dependence on flux. Fits for a variable emissivity 
index (M2-qv-Fe) and the l am ppost geometry (M2-LP-Fe), 
can be found in Tables |A7| and |A8| respectively. Note that 
a variable emissivity results in best-fit values consistent with 
q^3. 

As can be seen in Table 1^ the inclination is higher than in 
our other models, particulmy compared to M2b-LP and M3- 
q3. However, it is still consistent with being below 60° as 
required by the lack of eclipses and consistent with all M2 
models. See Section lAll for a more detailed discussion about 
the inclination. Furthermore, this is the only model where we 
find a significant truncation of the inner accretion disk around 
30 Tg. None of the other model parameters changes signifi¬ 
cantly. 

4. DISCUSSION 

We have presented simultaneous spectral fits to five 
Swfft/XRT and NuSTAR observations during the failed 2013 
outburst of GX 339-4 covering luminosities between ~0.9- 
6% LEdd- All observations show a very hard power-law, with 
a photon index F =» 1.6 and clear evidence for reflection. We 
have shown that standard models, consisting of a power-law 
continuum and an additive reflection model fail to reproduce 
the data within sensible physical parameters. The data, thanks 
to the very high S/N, clearly indicate that the input to the re¬ 
flection model needs to be different from the continuum emis¬ 
sion observed directly. Our best fit solutions typically require 
that the inner radius of the accretion disk extends close to the 
ISCO (Rin 100 rg). 

For com parison, we analyzed XMM-Newto n data pre¬ 
sented by ]Plant, O’Brien & Fender] (|2014|l taken dur¬ 
ing the same outburst (between our observations IV and 
V; XMM-Newton ObsIDs 0692341201, 0692341301, and 
0692341401). We used a similar annular extraction region 
as these authors with an inner radius of 11.5" and an outer ra¬ 


dius of 45" to excise the heavily pil ed-up core. Applying the 
best-fit relxill model presented by|Plant, O’Brien & Fender] 
(|2014|) and fitting all three observations simultaneously results 
in an acceptable fit with ;^f^gydof= 1.25/5079. When applying 
model M2-q3, and fixing the secondary power-law and the 
ionization to the best-fit NuSTAR values, we find a very sim¬ 
ilar quality of fit with ;^f^^^/dof= 1.26/5081 and similar struc¬ 
tures in the residuals. This similarity shows that the soft X- 
ray band-pass of XMM-Newton is not sufficient to constrain 
the complex accretion geometry, as the Compton hump is not 
covered and that at the same time our model is fully compati¬ 
ble with the XMM-Newton data. 

Besides the input power-law to the reflector, the biggest dif- 


ference between our models and the ones presented by Plant, 
O’Brien & Fender] (|2014)l is that they find evidence for the 
presence of a cold disk thermal component. Such a compo¬ 
nent isofteno^eryed in GX 339-4 in t he hard state (see, 
e.g, JMiller et al.|[2006[ |Reis et JI]|2008[ [Wilkinson & Utt-] 
ley||2009|l, but we do not find an improvement by adding it. 


Its absence is likely connected to the less sensitive XRT data 
compared to XMM-Newton as well as a degeneracy between 
the absorption column and the disk. This degeneracy is seen 


when describi ng the XMM-Newton data of Plant, O’Brien & 
|Fender| ( f2014| l with the M2-q3 model, where we find slightly 
lower values of the absorption column (Ah = (0.664+0.006) x 


10"^"^ cm ^ compared to « (0.74 + 0.01) x 10' 


, 22 , 


"). 


We note that the calculated reflection fractions are all below 
unity. In a static geometry, the reflection fraction corresponds 
approximately to the solid angle covered by the reflector in 
units of It:, i.e., a reflection fraction of 1 would correspond 
to an infinite slab illuminated from above. We measure val¬ 
ues below 0.5, which could indicate a truncation of the accre¬ 
tion disk at the inner radius. On the other hand, a outflowing 
corona would also result in a reduced reflection fraction, as 
the coronal radiation is beamed away from the accretion disk 
( |Beloborodov|1999] l. 

The photon index of the observed primary continuum seems 
to show a hysteresis effect, with observation V showing a sig¬ 
nificantly softer spectrum than observations I and II, despite 
being much fainter. This hysteresis could be connected to the 
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hyster esis observed in the radio/X-ray and N IR/X-ray corre¬ 
lation (|Corbel et al.|2013||Russell et al.|2007|l which might be 
related to different jet behavior between rising and decaying 
hard state observations. However, this hysteresis is typically 
observed in a full outburst, in which the hard states are sepa¬ 
rated by a soft or intermediate state, where we expect stronger 
changes in the accretion geometry. 

With our complex corona models (M2, M2b, and M3), we 
find an iron abundance of typically ~1.8 times solar, in line 


with previo us work and as expected in LMXBs (e.g. |Allured| 
et al.|201^. The fit statistically constrains the iron abundance 


extremely well, but the value strongly depends on the model 
assumptions. For example, with a variable emissivity index 
(models M2-qv and M3-qv) we find an iron abundance around 
2.5 times solar, while it drops to around 1.5 solar in all mod¬ 
els when a narrow Gaussian line is included (e.g., M2-q3-Fe). 
The absolute value and the error bars are therefore model de¬ 
pendent and subject to systematic uncertainties not accounted 
for in the tables. 


4.1. Accretion geometry 

As discussed in [Fabian et al.| ( |2014| l, due to light-bending, 
the regions of the corona closer to the black hole contribute 
more irradiating flux on the accretion disk than regions fur¬ 
ther away, while the opposite is true for the primary observed 
continuum. A change in the coronal parameters with height 
will therefore result in an observed continuum with a different 
spectral index than the one incident to the reflector. This effect 
is most relevant if the corona extends close to the black h ole 
where relativistic effects are strongest (jDauser et al.|2013|l. 

We can best approximate this geometry by using model 
M2b-LP consisting of two coronae at different heights above 
the black hole with different photon indices, as presented in 
Table H We find that the inner corona originates very close 
to the black hole, and thus is subject to strong light bending 
effects. These effects prevent most of the flux from that part 
of the corona from reaching the observer directly, and we only 
see the reflected part. 

We also find evidence for a complex shape of the Fe Kor 
line, which is clearly relativistically broadened with an ad¬ 
ditional narrow core, close to 6.4 keV. A detailed discus¬ 
sion of its shape is beyond the scope of this paper and will 
be presented in a forthcoming publication making use of 
Suzaku/XlS data with better spectral resolution (Tomsick et 
al., in prep.). The broadened component is constrained in our 
fits through the relconv smearing kernel. In the models with¬ 
out an extra component for the narrow core, the presence of 
a reflector within ^tlOOr^ is clearly required (see Table |^. 
Adding the narrow core moves the lower limit of the inner 
accretion disk radius out to about 20 rg. 

For the faintest observation (V) all models indicate the pos¬ 
sibility that the accretion disk is truncated. These values are 
only weakly constrained and strongly dependent on the as¬ 
sumed X-ray source geometry (see Tables|^and[^. It is clear, 
however, that during observation V the source spectrum was 
different from the first four observations (see also Figure 
with either a change in the inner accretion disk radius or the 
corona’s location and spectrum or a combination thereof. 


4.2. Inclination and spin 

All our models indicate an inclination in the range around 
50°. This value is significantly higher than i « 20-30° 
found when modeling other observations with similar re- 


flection models ([Miller et al.||2004t [Reis et al.||2008| |Plant,| 


[O’Brien & Fender|2014[ l. Kolehmainen & Done[ ( [20l6| l find 

that inclinations i > 45° give better fits in the high-soft state 
when fitting the disk continuum to measure the spin. Our re¬ 
sult seems to reconcile the continuum fits with the reflection 
fits with respect to the disk inclination. All of these previous 
data, however, were taken during a full outburst, i.e., an out¬ 
burst that followed the standard evolution through the high- 
soft state. The data presented here were taken during a failed 
outburst, which might have a different geometry. 

Our models do not const rain the spin . We use a - 0.92 
throughout, as measured by [Miller et al.[([2004|l. However, as 
the inner radius is only weakly constrained and typically of 
the order of 10 Cg, lower spin values are completely consis¬ 
tent with our results. To test this, we set the spin to a = 0 
for th e M2 -q3 model and obtained basically identical fits (see 
Table[A9ll. 


5. SUMMARY 

While the combined NuSTAR and Swift data provide one 
of the best data-sets on GX 339-4 in the low/hard state to 
date, we have shown that it is difficult to measure the inner 
truncation radius of the accretion disk precisely. The mea¬ 
sured value depends strongly on the assumed geometry and 
emissivity profile of the accretion disk. However, we find 
no evidence for a strongly truncated disk, i.e., with an in¬ 
ner radius > lOOr^. Furthermore, our spectral fits clearly 
show that the continuum spectrum incident to the reflector is 
significantly different from the observed primary continuum. 
The spectrum reflected by the accretion disk is significantly 
harder, which is necessary to explain the relative strength of 
the Compton hump to the Fe Ka line. A lamppost geome¬ 
try with changing spectral hardness as a function of coronal 
height seems to describe the observed spectra well, but can 
only be regarded as a crude approximation to the true physi¬ 
cal geometry. 

We would like to stress again that the data were taken dur¬ 
ing a failed outburst, during which the source did not switch 
into the soft state. It is currently unknown what the difference 
between failed and standard outbursts is, and how a transi¬ 
tion to the soft state is triggered. Continued monitoring of 
GX 339-4 and similar black hole transients is necessary to 
answer these questions and study if we can measure signif¬ 
icantly different accretion geometries in these two types of 
outbursts. 
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TABLE Al 

Best-fit parameters for an emissivity index q = 2> and spin a = 0.92 (M1-q3). 


Parameter 

I 


II 


Ill 


IV 


V 

Ah [1022cm-2] 

A 7O+J+0.020 

U./83-0 016 

— 


— 


— 


— 


Fe/solar 
i [deg] 

J.UU-Q ,2 

3Q 0+0 / 

-1.2_ 

— 


— 


— 


— 



“01620^™^ 
1.562 ±0.004 
2.41+«:[3)x10-^ 


^cont^ 

r 

Arefl" 

Rm ['"si 
CCxRT 


0.0607 ± 0.0007 
1.540 ±0.007 

0.97+“ 21) X 10-5 


2.04!“ >2) X 102 

4 6 +°“ 

+^-0.6 

0 962+“ “27 
0 078 


0.0999 ± 0.0008 
1.544 ±0.005 
1.67+“ |“)x 10-5 


2.11 


+0.08 


0.06 

A 9+0.6 

-0.5 

n.oi 5 


X 10^ 


2.32+“:‘>)x102 

4-0!oi 


^•^^'^^-O.OOM 

1.587 ±0.004 
3.57+“;W) X 10-5 


2.20!“:“2) X 102 
T SQ+^4o 
1.138_„„P 


0.0391 ± 0.0004 
1.573 ±0.006 
(4.6 ± 0.4) X 10-“ 
(2.09!“;“*) X 102 

A ^+0.8 

-0.5 

i.UUi-0_^14 


R 

().48 

0.54 

039 

(154 

0.42 


1.63 

2.72 

4.26 

5.49 

0.92 

^2/d.O.f. 

+?ed 

2391.96/1703 

1.405 






‘ in ph s ^ cm ^ 

’ Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of 10 Mq 


TABLE A2 




Best-fit parameters for a free emissivity index 

(MI-qv). 


Parameter 

I 

II 

III 

IV 

V 

Ah [1022cm-2] 

0.794 ±0.019 

_ 

_ 

_ 

_ 

Fe/solar 

6.5 ± 0.4 

— 

— 

— 

— 

i [deg] 

43.0/1-5 

— 

— 

— 

— 

Acont^ 

0.0611+“-“““' 

0.1008 ±0.0008 

0.1639 ±0.0010 

0.2289 ± 0.0012 

0.0390 ± 0.0004 

r 

1 545+a8-r5 

1.549 ±0.005 

1.569 ±0.004 

1.593 ±0.004 

1.572 ±0.006 

Arefl" 

(i.07!“;55)x10-5 

(1.81 ±0.10) X 10-5 

(2.63 ± 0.15) X 10-5 

(3.90 ± 0.18) X 10-5 

(0.47 ± 0.04) X 10-5 


(2.04/“;>j)xl02 

(2.14/“;“«)x102 

(2.35/“:>“)x102 

(2.21/“:““) X 102 

(2.12/“:“5) X 102 

Rin ['■j] 

+3/p 

5 0+2-“ 

+ ■ -1.0 

4.4/“-“ 

+2/!ii 

4.0 ± 0.4 

4.0/“:“ 

q 

6 2+2-5 

5 4+* -“ 

•7- -0.7 

4 8+“ “ 

1.7 ±0.5 

CCxRT 

0.964 ± 0.028 

0.935 ±0.015 

1.161 ±0.014 

1.138 ±0.012 

1.065 ±0.014 

R 

0.52 

0.58 

0.63 

0.68 

0.43 

%Ledd'’ 

1.62 

2.71 

4.25 

5.48 

0.93 

^2/d.O.f. 

2214.67/1698 





+L 

1.304 






® in ph s ^ cm ^ 

^ Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of 10 M© 


TABLE A3 

Best-fit parameters for the lamppost geometry (Ml-LP). 


Parameter 

I 

II 

III 

IV 

V 

Ah 

Fe/solar 

([deg] 

0.786/“||« 

c 9'7+0.37 
-0.29 

40 

E 

E 

E 

E 

Acont'* 

r 

A(corona)[rj] 

Rin[fj] 

CCxis 

““Q“QgQ9+(IDDD7 

^•^^^^-0.0005 

i.j‘+z,-ooo5 

(9.9544939/“:—)x 10-“ 

( 2 . 04 /“:“) X 102 

2 1+2-' 

46/8:8 

-0.5 

027 

““Q'^^A2+tn)D08 

-0.0007 

^•'^^^-0.004 

(1.71 ±0.10) X 10-5 
( 2 . 11 /“:“ 2 )x 102 

9 1 + 0.8 
- 0.0 

4.3 ± 0.4 

0.936 ±0.015 

0 1626+“-““'“ 

1.565 ±0.004 

(2.47/“:|5)x10-5 

(2.3243/“:“»g)xl02 

9 1 + 0.6 

4'i6+«2“ 

^•^'^-0.32 

0.2275 ± 0.0012 
1.5893/“:““55 
(3.65/“:|2)x10-5 
(2.2078/“:“55“) X 102 

2 l+“-5 

4 5+8-25 

1 1 97 + 0.012 

0.0391 ± 0.0004 

1.574 ±0.006 
(4.7 ± 0.4) X 10-“ 
(2.09/“:“2Jx102 

2 1 + 2-2 

4 i/8:8 

^•^-0.5 

1.061 ±0.014 

R 

"1149 

“1)35 

“ira) 

“1155 

“1143 

LJUdd X 100 “ 

1.63 

2.71 

4.25 

5.51 

0.92 

4-2/d.o.f. 

2330.59/1698 





xid 

1.373 






^ in ph s ^ cm ^ 


Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of 10 M© 
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TABLE A4 

Best-fit parameters for a free emissivity parameter and untied photon indices (M2-qv). 


Parameter 

I 


II 


Ill 

IV 


V 

Ah [lO^^cm-^] 

0.826!®;“® 

— 


— 


— 

— 


Fe/solar 
i [deg] 

9 09+0.25 

44 3+2-^ 

— 


— 


— 

— 



A A 

^cont 

Tpower 

Trefi 

Arefl" 

f 

[rs] 

CCxRT 


0.0655 ± 0.0010 

, =74+0.009 
-0.008 
1.33 ±0.05 

l.OC^^^lx 10-5 

2.09!® >5) X 10^ 
4 0+14.3 
’- 2.8 

0.949+^028 


^•^^^^-0.0012 
1.584 ±0.007 

1 '2'37+0.031 

i.jj/_0 028 

1.73!® |2jx 10-5 

2-2lti?)xl02 

2!o 

0.917 ±0.015 


1.605 ±0.006 
1.364 ±0.020 
2.33!« [^) X 10-5 
2.53!®:|^)x102 
2+6 

l.B4+°®!! 


^ -0.0027 

1.633 ±0.005 

1 oqc+0.014 
i.j5'.7_q018 

3.40!«;>nxl0-5 
2.40!®;>« X 102 
2-1!^:^ 

1.108+°°12 


0.0412 ± 0.0005 
1 596+®'®®’ 

i.j+U-0 006 
1.37 ±0.05 
(4.2 ± 0.4) X 10-® 
(2.25!«;‘2)x102 
2+81 

1.052+;!-014 


AL 

0.25 ± 0.05 

0.247 ± 0.030 

0.241 ± 0.021 

0.239 ± 0.016 

0.23 ± 0.05 

R 

0.60 

0.68 

0.73 

0.79 

0.48 


1.63 

2.72 

4.30 

5.55 

0.93 

4^9d.o.f. 

+?ed 

2077.77/1693 

1.227 






' in ph s ^ cm ^ 

’ Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of 10 Mq 


TABLE A5 

Best-fit parameters for a free emissivity parameter and two power-law continua (M3-qv). 


Parameter 

I 

II 

Ill 

IV 


V 

Ah [10“cm-2] 
Fe/solar 

0.884!®;- 

9 C9+O.26 
^•‘^^-0.20 

— 

— 

— 

— 


i [deg] 

36!l 

— 

— 

— 

— 



< 0.6 X 10"^ 


■^cont, 1 
rcont,l and refl 
Arefl" 

^in [fg] 

q 

^cont,2^ 

rcont,2 

CCxRT 


0 

-QQ31 

1 47J71+0.022 

0.66!0 ‘o X 10-5 


^•^^'^-0.008 
1 425+0-01 / 

1.23+Slf)x 10"^ 


2.38+®;g) X 10^ 
0+44 
°-6 
1 4+0 9 
4-0.7 

0035+®®04 


2.42!®-}^) X 10^ 


1 . 
0.931 


0.005 

+0.07 


911 0+12.295 
^••^^'^-0.004 
171 +0.42 
-0.23 

^•^uj_o.oio 
1 R8+^'^^ 

0:904?- 


1 422+()-0iu 

i.'t^j-0 0i7 

(l.64+[j;3)xl0-^ 

(2.95 ± 0.20) xiO^ 

9+10 

^-0 

1 cg+0.23 
^•‘^^-0.24 
0.116 ±0.008 
1 SSS+O-059 

1.0J0_QQ3Q 

1.123+Sa1^ 


0 091+^^»^4 

-QQ19 
1 441+0.009 
^•^^^-0.010 
2.51+®ll)x 10-5 


2.77!«;25)x102 

9 191 +4.017 
^•■^^•^-0.012 
I 90+019 

-0.12 

0.170!«;«|^ 

i-834!“;“« 

1.102 ± 0.010 


1 IQI+0.029 
A-’^-^-0.040 
(3 59+0.30) 

IO.OV_O.34) 

(2 50‘^^-^^ I 
IZ.jU q 13 I 

- 0 . 15 ) 

I 4 + 8.7 
^ -0.9 

0 0024 

^•‘^^^-0.007 

1.067+®®12 


X 10-® 
X 10^ 
X 10^ 


AL 

0.45 ± 0.06 

0.45 ±0.12 

0.43 ± 0.05 

0.39 ± 0.04 

0.21 ±0.04 

Rl 

0.48 

0.57 

0.67 

0.84 

INF 

Ri+2 

0.32 

0.36 

0.38 

0.43 

0.48 

%iedd® 

2.11 

3.57 

5.76 

7.52 

0.94 

4-9d.o.f. 

2043.38/1688 





xii 

1.211 






^ in ph s * cm ^ 

® Luminosity calculated between 0.1-300keV assuming a distance of 8kpc and a black hole mass of lOM© 
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TABLE A6 

Best-fit parameters for the lamppost geometry and two power-law continua (M3-LP). 


Parameter 


I 


III 


IV 


Nu [10^2 cm-2] 

Fe/solar 
//(corona) [r^] 
'■[deg] 


0.908!«p 

44+23 

3ii 


+0.004 

-^.p06 


^+0.036 
-0.020 
1 a c+0.06 
^••^'^-0.08 
1 433+0.046 


^0T20gSl 

1.486t“i 


“TrT43^RnK)9 

1 495+0.021 
i.‘t^O_0 084 

1 431+0.013 

-0.014 

(2.30 ± 0.24) X 10^5 

(2.87^° 34) X 102 
2-2005!i'j,“3 

n 190+0.072 

U-1^0-0 028 

2 OO+® '7 

-0.14 

1.098 ±0.013 


l-21ao 

1 41 + 0.04 
^•^^-0.06 

) X 10"^ 


A 

2 ^cont 

Tpower 

Trefl 

Arefl" 

Rin [rg] 

A + 1 ^ 

2 ^cont ,2 

rpower,2 

CCxis 


0.059 

1 72+^-^ 

^•'^-0.04 
1 4 C+ 0 .O 8 

(0-55t‘° X 10-3 


0.035+ 


2-49tg)xl02 

(l.4+[S)x 102 

0 011+0 028 
+'-+'^^_q.006 

1 94 + 0.18 
0933+8^25 

u.yjj -0(^28 


-0.024 

1.06t|«)xl0-3 


2.51+^-23)x 10^ 

99+24 
^^-28 
n n8T+o oo6 

^•0o3_qq^q 

1 09+0.09 
^•°"^- 0 .p 8 

u.vuH -0017 


3.03+|!;g) X 10^ 

40+26 

^^-38 

u.u/y-0 012 

2 10 ^*^ 

-0.14 

1.1I8±0.018 


W- -o!4j 

I 9 ^7+0.18 
-0.15 




X 102 
X 102 


0.0408!°;™^ 
1 fiT+o^ 
1066+®‘3 

i.UOO -0(^l4 


AL 

0.49 ± 0.21 

-0.48 ±0.11 

-0.61 ±0.15 

-0.51 ±0.17 

-0.4 ± 1.0 

R 

0.92 

0.68 

0.63 

0.72 

1.36 

L/Ldd X 100® 

1.74 

2.91 

5.03 

5.85 

0.96 

^2/d.O.f. 

xid 

2058.22/1686 

1.221 






‘ in ph s ^ cm ^ 

’ Luminosity calculated between 0.1-300keV assuming a distance of 8kpc and a black hole mass of lOM© 


TABLE A7 

Best-fit parameters for a free emissivity parameter, untied photon indices and an additional narrow Gaussian line at 

6.4 keV (M2-qv-Fe). 


Parameter 

I 


11 


III 


IV 


V 

Ah [1022cm-2] 
Fe/solar 
! [deg] 

A occ+0.020 
U.825_gQjg 

^•‘^^-0.09 

60+J« 

= 


E 


E 


E 



Acont 

rcont 

Trefl 

Arefl' 


0.0667 ± 0.0010 0.1116 ± 0.0014 


1 SST+^-^^^ 
i.joO-0 005 

1.24 ±0.06 

(1 07+°®®’! 

( + ■+^'-0.121 

2.10tlh 

f0.48t32) 


X 10-3 
X 102 
X 102 


1 5Q5<;+00073 

1.27 ±0.04 

I X 10-3 

|X 102 

lx 102 


Cl 7g+0.13 

G-'^-0.12, 

12 24+0-09 


^•^^^*^-0.0019 
i.ui /-7 _qqq29 
1 979+0.028 
^•"^'^-0.026 
(n 72Q+0.14 


q 3 


ApeKo-^ 

GCxrt 


+6.3 
2.1 

, Q9+0.28 

• -0.32, 

0.950 ± 0.028 


(0.1 


X 10- 


2-4!l;2 

(0.7 ± 0.5) X lO-”* 
0.911 ±0.015 


2.59+«;i^ 

lo+S 


X 10-3 

X 102 
X 102 


^•^-^^^-0.0024 
1 6459 + 0-0048 

■3,"+d.8fr 

1.010 0017 
(3.44 ± 0.18) X 10- 


X 10+ 


1.124 ± 


I X 10-^ 

0.014 


(2-44l!’;i“) 

44+32 

3.2+f^ 

(2.1 ± 0.6) X 10-"^ 
1.099 ±0.012 


0.0418 ±0.0006 
1 603+®™* 

l.UUO pQQ, 

1.24 ±0.06 
(4.4 ± 0.5) X 10-3 
(2.32+0:|3)xl02 
(0.8 ± 0.5) X 102 
10 +® 

( 0 . 54 + 0 + 2 ) X 10 -* 

1.053 ±0.014 


AL 

0.35 ± 0.06 

0.33 ± 0.04 

0.345 ± 0.027 

0.330 ±0.018 

0.36 ± 0.06 

R 

0.66 

0.74 

0.81 

0.86 

0.53 


1.64 

2.74 

4.37 

5.62 

0.94 

X^/d.o.i. 

2005.72/1688 

1.188 






^ in ph s * cm ^ 

Luminosity calculated between 0.1-300 keV assuming a distance of 8 kpc and a black hole mass of 10 Mq 
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TABLE A8 

Best-fit parameters for the lamppost geometry, untied photon indices and an additional narrow Gaussian line at 6.4 keV 

(M2-LP-Fe). 


Parameter 

I 

II 

III 

IV 

V 

Ah 
F e/solar 
A(corona) [r^] 
i [deg] 

0.853!°;“^ 

1 ^C+h.lO 

— 

— 

— 

— 

^tcom 

Tpower 

Lefl 

Rin [rg] 
AFeKa^ 
CCxRT 

0.0667 ±0.0011 
1.583 ±0.009 

1.24 ±0.06 
fl.08+|‘';1x 10-® 
(2.10+0;‘3 )x102 
39+22 

(7.8;2.3^ X 10-® 

0.949 ± 0.028 

0.1114 ±0.0014 

1 5Q44+0.0068 
-0.0011 

1.27 ±0.04 

fl.75+‘; |2i X 10-® 
2.24 +<;;o9) X 102 
( 7 +®) X 10 ' 

(7 ± 4) X 10-® 

A qi 9+0.014 

0.1823 ±0.0018 

1 ^io+0.006 

'■"'»-0.0Q4 

1 97Q+0.026 

^•"^'^-0.027 

(2.40 ± 0.16) X 10-® 
(2.58!«;i^)xl02 

( 1 . 7 +ohx 10-^ 

1.124 ±0.014 

0.2556 ± 0.0023 

1.645 ±0.005 

1 31 Q+ 0.017 

(3.45 ± 0.17) X IO-® 

(2-44t1)xl02 

40t\l 

(2.1 ±0.5) X 10-^ 
1.099 ±0.012 

0.0418 ±0.0006 

1 603+®™ 
i.uuj PQ07 

1.25 ±0.06 
(4.4 ± 0.5) X 10-® 

(2-3i1:lo)xio" 

45 +g 

(5.1+|'2)x 10-® 
1.052 ±0.014 

AT 

0.34 ± 0.06 

0.33 ± 0.04 

0.338 ± 0.027 

0.327 ± 0.017 

0.36 ± 0.06 

R 

0.66 

0.73 

0.80 

0.85 

0.53 

Eqw (EeKa) [eV] 

21.22 

11.59 

17.33 

16.17 

23.18 

L,/L,dd X 100® 

1.64 

2.74 

4.36 

5.62 

0.94 

X^/d-o.f. 

2008.80/1692 





aL 

1.187 






^ in ph s * cm ^ 

^ Luminosity calculated between 0.1-300keV assuming a distance of 8kpc and a black hole mass of lOM© 


TABLE A9 

Best-fit parameters for emissivity index q = 3, spin <3 = 0, and untied photon indices (M2-q3-a0). 


Parameter 

I 


II 


Ill 


IV 


V 

Ah [1022cm-2] 
Fe/solar 
i [deg] 

0.868 ± 0.020 

1 7T+h.09 

A-'-’-0.08 

47+;2 

— 


— 


— 


— 



A 

^cont 

Tpower 

Trefl 

■4refl 


L'gJ 

GCxrt 


^•^'^'^^-0.0013 

1.585 ±0.009 

1 2Q+0.07 

(9.6!{ 3) X 10-® 
)x 10^ 


0.1112 ±0.0014 
1.594 ±0.007 

1 Q 19+0.025 
^•■^^^-0.033 


1 . 66 + 




S+0.12 
-0.11 
['2 28+0 09 


|x 10^ 
0.945 ± 0.028 




-0.08 
|x 10' 
0.915 ±0.015 


X 10-® 
X 10^ 


(91) 


0.1812 ±0.0019 
1.617 ±0.006 
1 099+0.017 
i.j.:?o_00i6 

(2.25 ± 0.14) X 10- 

2.65+1; ‘^) X 10^ 

1.3!J:«)x102 
1.129 ±0.014 


0.2538 ± 0.0023 
1.643 ±0.005 

1 9C-7+O.OI6 
i.OJ/_00i5 

3.23t|hxlO-® 


2-52l:!‘)xl02 
65+g 

1.105 ±0.012 


0.0420 ± 0.0006 
1.608 ±0.007 
1.34 ±0.04 
(4.2 ± 0.4) X 10-® 
2-34l:!?)xl02 

2.3;° h X 10^ 

1.052 ±0.014 


AL 

0.30 ± 0.06 

0.282 ± 0.029 

0.284 ±0.018 

0.285 ±0.016 

0.26 ± 0.04 

R 

0.64 

0.71 

0.78 

0.83 

0.52 

%iedd® 

1.65 

2.74 

4.35 

5.63 

0.93 

T+d.o.f. 

xii 

2123.79/1698 

1.251 






^ in ph s ^ cm ^ 

Luminosity calculated between 0.1-300 keV assuming a distance of 8kpc and a black hole mass of lOM© 



